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ABSTRACT 

Context. In the project Galactic Cold Cores we are carrying out Herschel photometric observations of cold regions of the interstellar 
clouds as previously identified with the Planck satellite. The aim of the project is to derive the physical properties of the population 
of cold clumps and to study its connection to ongoing and future star formation. 

Aims. We examine the cloud structure around the Planck detections in 71 fields observed with the Herschel SPIRE instrument by the 
summer of 2011. We wish to determine the general physical characteristics of the fields and to examine the morphology of the clouds 
where the cold high column density clumps are found. 

Methods. Using the Herschel SPIRE data, we derive colour temperature and column density maps of the fields. Together with 
ancillary data, we examine the infrared spectral energy distributions of the main clumps. The clouds are categorised according to their 
large scale morphology. With the help of recently released WISE satellite data, we look for signs of enhanced mid-infrared scattering 
('coreshine'), an indication of growth of the dust grains, and have a first look at the star formation activity associated with the cold 
clumps. 

Results. The mapped clouds have distances ranging from ~100pc to several kiloparsecs and cover a range of sizes and masses 
from cores of less than 10 Mq to clouds with masses in excess of 10000 Mq. Most fields contain some filamentary structures and in 
about half of the cases a filament or a few filaments dominate the morphology. In one case out of ten, the clouds show a cometary 
shape or have sharp boundaries indicative of compression by an external force. The width of the filaments is typically ~ 0.2-0.3 pc. 
However, there is significant variation from 0.1 pc to 1 pc and the estimates are sensitive to the methods used and the very definition 
of a filament. Enhanced mid-infrared scattering, coreshine, was detected in four clouds with six additional tentative detections. The 
cloud LDN 183 is included in our sample and remains the best example of this phenomenon. About half of the fields are associated 
with active star formation as indicated by the presence of mid-infrared point sources. The mid-infrared sources often coincide with 
structures whose sub-millimetre spectra are still dominated by the cold dust. 
Conclusions. 

Key words. ISM: clouds - Infrared: ISM - Submillimeter: ISM - dust, extinction - Stars: formation ~ Stars: protostars 



1. Introduction 

The main phases of the star formation process are largely un- 
derstood, starting from molecular clouds a nd progressing via 
dense cores down to protostellar collapse (iMcKee & Ostriken 



* Planck {http://www.esa.int/PlmcQ is a project of the European 
Space Agency - ESA - with instruments provided by two scientific 
consortia funded by ESA member states (in particular the lead coun- 
tries: France and Italy) with contributions from NASA (USA), and tele- 
scope reflectors provided in a collaboration between ESA and a scien- 
tific Consortium led and funded by Denmark. 

** Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



I2OO7I) . This is thanks to the detailed observations of the near- 
est low mass and intermediate mass star formation regions and, 
on the other hand, sophisticated numerical modelling. However, 
the formation of each star is an individual process and, for a 
global view, extensive surveys of different environments and of 
the different phases of the star formation process are needed. 

One major question is how the properties of star forma- 
tion depend on the initial conditions within the cold molec- 
ular clouds. Direct links should exist between the physics of 
the prestellar clumps and the star formation efficiency, the 
mode of star formation (clustered v s. isolated), and the masses 
of the bom stars (j Elmegreen 20 ij; IPadoan & Nordlundll201 U 
iNguven Luong et al..,201 1.) . The general time scales of star for- 
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mation have been extensively studied. Eventually, one has to 
progress to the examination of the finer details, how the pro- 
cesses vary between different environments and, in each case, 
what is the interplay between turbulence, magnetic fields, kine- 
matics, and gravity. Similarly, our understanding of the role of 
external triggering is incomplete. How does the influence of su- 
pernovae and stellar winds (and of galactic density waves, col- 
liding flows of HI gas, etc.) vary between different regions of the 
Milky Way and is this reflected in the timescales and the stellar 
initial mass function? 

The Planck satellite (iTauber et al.ll20ldh has made a new ap- 
proach to the study of the earliest stages of star formation possi- 
ble. By mapping the whole sky at several sub-millimetre wave- 
lengths, with high sensitivity and small beam size (below 5' at 
the highest frequencies), Planck is providing data for a global 
census of the coldest component of the interstellar medium. 
The selection of the cold (Tdust <14K) and 'compact' (close to 
beam size) objects has led to a list of more than 10000 objects. 
Because of the limited reso lution, this Cold Clump Catalogue of 
Planck Objects (C3PO, see' Planck Collaboration et al.ll201 id) is 
not dominated by cores (pre-stellar or only starless cores at sub- 
parsec scales) but by ~1 pc sized clumps and even larger cloud 
structures extending to even tens of pc in size. However, the 
low temperatures of the objects (below 14 K) are possible only 
for the denser, less evolved regions that are well shielded from 
the interstellar radiation field. Therefore, the objects detected by 
Planck are likely to contain one or several cores, many of which 
will be pre-stellar or in the early stages of protostellar evolu- 
tion. This Planck survey constitutes the first unbiased census (in 
terms of sky coverage) of possible future star forming sites and 
provides a good starting point for global studies addressing the 
pre-stellar phase of cloud evolution. 

Within the Herschel Open Time Key Programme Galactic 
Cold Cores, we are mapping selected Planck C3PO objects 
with the Herschel PACS and SPIRE instruments (100-500 Aim). 
The higher spatial resolution of Herschel ( Poglitsch et aL. 2010t 
iGriffin et al.l 1201 Ol) makes it possible to examine the structure 
of the sources which gave rise to the Planck detections, often 
resolving the individual cores. The inclusion of shorter wave- 
lengths (down to lOOjum) help to determine the physical char- 
acteristics of the sources and their environment, and to investi- 
gate the properties of the interstellar dust grains. Our Herschel 
survey will eventually cover some 120 fields between 12' and 
one degree in size and will altogether cover approximately 
350 individual Planck detections of cold clumps. Preliminary 
results of th is follow-up have been presented in lJuvela et alJ 
(1201 Ol 1201 lb (papers I and II, respectively) from the Herschel 
Science Demonstration Phase observations (three fields there 
called PCC288, PCC249, and PCC550) and on a sam ple of ten 
Planck sources in lPlanck Collaboration et al] ( 1201 Ibl) . 

In this paper we present results for the first 7 1 fields that were 
observed with the SPIRE instrumenfl (250, 350, and 500 jum) by 
the summer 20 11. We concentrate on the large scale structure of 
the clouds and the general characteristics of the main clumps. 
We compare the sub-milUmetre emission with other available 

' SPIRE has been developed by a consortium of institutes led by 
Cardiff University (UK) and including Univ. Lethbridge (Canada); 
NAOC (China); CEA, LAM (France); IFSI, Univ. Padua (Italy); lAC 
(Spain); Stockholm Observatory (Sweden); Imperial College London, 
RAL, UCL-MSSL, UKATC, Univ. Sussex (UK); and Caltech, JPL, 
NHSC, Univ. Colorado (USA). This development has been supported 
by national funding agencies: CSA (Canada); NAOC (China); CEA, 
CNES, CNRS (France); ASI (Italy); MCINN (Spain); SNSB (Sweden); 
STFC (UK); and NASA (USA) 



far-infrared data, especially the AKARI (iMurakami et al.ll2007l) 
satellite far-infrared maps, and use Wide-field Infrared Survey 
(WISE Wright et al. 2010) satelUte observations (3. 6-22 //m) to 
look for regions where enhanced mid-infrared scattering could 
indicate an increase in the size of the dust grains. This so- 
called coreshine phen omenon was first dete cted with Spitzer 
data of LDN183 (see Stein acker et al.l 1201 0'). a cloud also in- 
cluded in the present sample. The point sources detected in the 
mid-infrared WISE data also serve as an indicator of the ongoing 
star formation. 

The structure of the paper is the following. The observa- 
tions are described in Sect. |2] The main results are presented 
in Sect. [3] starting with the analysis of the large scale properties 
of the clouds. These include the estimation of the colour tem- 
peratures (Sect. 13.1.1b . the column densities, and cloud masses 
(Sect. IlL^ . all based on the SPIRE observations. In Sect. [TO] 
we describe a categorisation of the target fields that is based on 
the main cloud morphology. In Sect. 13.21 we characterise the 
fields further with the help of infrared data. We examine the 
properties of the selected clumps, present their spectral energy 
distributions and masses (Sect. l3.2TTI ) and then look for coreshine 
in the mid-infrared data (Sect. [JX2b . The results concerning the 
general nature of the fields and the connection to the star forma- 
tion are discussed in Sect.|4] The final conclusions are presented 
in Sect. |5] The online appendices contain additional figures, in- 
cluding the surface brightness and the column density maps and 
the SEDs of selected clumps. 

A detailed analysis of all compact sub-millimetre clumps (in- 
cluding the core mass spectra) and of the young stellar objects 
associated with the cold clumps will be presented in future pa- 
pers. An in-depth study of dust emission properties (opacity and 
emissivity spectral indices) will also be presented later, together 
with the data from observations with the Herschel PACS instru- 
ment. 

2. Observations 

2.1. Target selection 

The Planck satellite is performing surv eys of the full sky a t nine 
wavelengths between 350yum and 1 cm dTauber et al.l2010 ). The 
analysis of the first two sky surveys led to the detection of more 
than 10000 compact sources of cold dust emission as described 
in Planck Collaboration et al. (201 Ic). The detections are based 
on the sub-millimetre cold dust signature that becomes visible 
when the warm extended emission is subtracted using the IRAS 
lOOjum maps as its spatial template a nd the average spe ctrum 
of the region as its spectral template (Mo ntier et al.ll2010 ). The 
detection procedure limits the maximum size of the detected 
clumps to ~ 12'. Together with the IRAS 100 /vm data, the Planck 
measurements have been used to estimate the colour tempera- 
tures of the sources. The source dis tances are e stimated using 
several methods (see P lanck Collaboration et alj|201 Id, for de- 
tails) that include association with known molecular clouds com- 
plexes or individual sources and the reddening of background 
stars at optical or near-infrared wavelengths ([Marshall et al.l 
2009; Mc Gehee 2012 in preparation). Altogether distance esti- 
mates exist for approximately one third of the C3P0 catalogue. 

The pre-selection of targets for the Herschel observations 
was based on a binning of Planck cold clumps with the fol- 
lowing parameter boundaries: Z= 0, 60, 120, and 180 degrees, 
\b\ =1, 5, 10, and 90 degrees, rdust=6, 9, 11, and 14K, and 
M =0, 0.01, 2.0, 500, 10^ M0. Here Tdust is the clump tem- 
perature obtained after the subtraction of the warm emission 
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d (pc) 

Fig. 2. The distance distribution of the fields. The bin at negative 
values corresponds to the sources without a distance estimate. 



component (see lPlanck Collaboration etal ]l2011ch. The binning 
was used to ensure a full coverage of the respective parameter 
ranges while at the same time weighting the sampling towards 
sources at high latitudes and with extreme values of the mass. 
The lowest mass bin (M-Q) was reserved for sources without 
distance estimates and, therefore, without any available mass 

o 

estimates. The Galactic latitudes \b\ < 1 were excluded be- 
cause those region s will be covered by the Herschel key pro- 
gramme Hi-GAL (iMolinari et alJl20Toh . Similarly, the regions 
covered by the oth er key progra mmes like the Gou ld Belt survey 
dAndre et alj2010l) or HOBYS dMotte et alj2010l) were avoided. 

The above procedure results in 108 bins, from which one 
target per bin was selected for Herschel follow up observations. 
Within each bin, the sources were inspected visually, comparing 
the Planck maps with IRAS and AKARI dust emission maps, 
extinction maps calculated using 2MASS star s, and CO emis- 
sion maps (Dame et al.l200lblFukui et al.ll999l) . when available. 
These data were used to confirm the reliability of the original de- 
tection. Some preference was given to fields containing several 
Planck detected clumps. The final selection of 108 fields cov- 
ers some 350 cold clumps from the C3PO catalogue. Out of the 
full list of 108 fields, in this paper we use the SPIRE observa- 
tions of the 71 fields for which the observations were completed 
by the Herschel observational day 721. In this sample, distance 
estimates exist for 75% of the fields (53 out of 71), partly be- 
cause a literature search has resulted in a few additional distance 
estimates. Figure [1] shows the positions of the sources on the 
sky and Fig. |2]the distribution of their distances. For the fields 
with distance estimates, the median value is 450 pc. When dis- 
tances were estimated from the stellar reddening ("Marsh all et al.l 
[2009; Mc Gehee 2012 in preparation), every effort was made to 
identify the extinction features most likely to be associated with 
the main clouds. However, the fields may contain clouds also 
at other distances, especially at the low Galactic latitudes. The 
distance distribution of all Planck clumps is smooth over the 
shown range of distances (Planck Collaboration et al. 2011c). 
The small number statistics and the uncertainties of the distance 
estimates (e.g., the use of different methods in different distance 
ranges) may contribute to the dips seen in Fig.|2]at ~600pc and 
~1400pc. The distribution is also affected by the exclusion of the 
areas covered by the other Herschel key programmes. In partic- 
ular, the avoidance of the Gould Belt clouds reduces the number 
of objects close to 500 pc. 



2.2. Herschel observations 

In this paper we use data from the first 71 fields observed 
with the Herschel SPIRE instrument. The observations con- 
sist of 250, 350, and 500 //m maps with the map sizes rang- 
ing from 37 to 77 arc minutes. The observations of the first 
three sources were done as part of the Science Demonstration 
Phase in November and December 2009 and the latest were 
performed at the beginning of May 2011. The fields are listed 
in Table [1] The Herschel observations were reduced with the 
Herschel Interactive Processing Environment HIPE v.7.0. using 
the official pipelinfl The resulting SPIRE maps are the product 
of direct projection onto the sky and averaging of the time or- 
dered data (the HIPE naive map making routine) In order of 
wavelength, the resolution of the SPIRE maps is 18", 25", and 
37", respectively. 

We rely on the Herschel calibration of the data. The ac- 
curacy of the absolute calibration of the SPIRE observations 
is expected to be better than 7'%|3 However, in order to de- 
termine the absolute zero point of the intensity scale, we car- 
ried out a comparison with Planck satellite observations that 
were complemented with the IRIS version of the IRAS lOOyum 
data (Miville-Deschenes & Lagache 2005). The Planck and IRIS 
measurements were interpolated to the Herschel wavelengths us- 
ing fitted modified blackbody curves, Bv(rdust)v^, with a fixed 
value of the spectral index, /3 = 2.0. The linear correlations be- 
tween Herschel and the reference data were extrapolated to zero 
Planck (+IRIS) surface brightness to determine the offsets for 
the Herschel maps. The uncertainties of the offsets were ob- 
tained from the formal errors of these fits. The errors are typi- 
cally ~1 MJy sr"' at 250 yum and, corresponding to the decreas- 
ing surface brightness values at the longer wavelengths. The de- 
rived intensity zero point is independent of the Planck calibra- 
tion and any multiplicative errors included in the comparison. 
Similarly, the offset determination depends very little on the as- 
sumed value of j6. Colour corrections can have an effect on the 
quality of the linear correlation between Planck and Herschel 
data. The data were first colour corrected assuming a constant 
colour temperature of 15 K. The offsets were later re-evaluated 
using the actual dust colour temperature estimates (Sect. [mi l. 

2.3. Other infrared data 

In the infrared range, we use data from the AKARI and WISE 
satellites. 

From the AKARI survey (iMurak ami et al.ll2007h we use ob- 
servations from the all-sky survey made with the FIS instrument. 
These include observations in the narrow band filters N60 and 
N160 (central wavelengths 65yum and 160 yum) and in the wide 
band filters WideS and WideL (central wavelengths 90 fim and 
140 jum). The spatial resolution of the data ranges from 37" at 
65 /im to 61" at 160jum. The accuracy of the calibration is as- 
sumed to be ~30%. AKARI data are available for all the fields. 

The WISE satellite (Wright et al. 2010) has four bands cen- 
tred at 3.4, 4.6, 12.0, and 22.0yum with spatial resolution rang- 

- HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 

^ The reduced data will become available through the ESA web 
site http://herschel.esac.esa.int/UserReducedData.shtml Further infor- 
mation is available on the home page of the Cold Cores project 
https://wiki. helsinki.fi/display/PlanckHerschel/The+Cold+Cores. 

SPIRE Observer's manual, 
\http://herschel.esac.esa.int/Documentation.shtml\ 
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Fig. 1. The locations of the observed fields. The background is the IRAS 100 //m map and the circles denote the positions of the 
fields. The arrows drawn on each circle indicate the distances, starting with Opc for the upright direction, one clockwise rotation 
corresponding to 2kpc. The six sources with distances larger than 1.5 kpc are drawn with smaller symbols. The sources without 
reliable distance estimates are marked with crosses. 



ing from 6.1" at the shortest wavelength to 12" at 22yum. The 
first public release of WISE satellite data took place in April 
2011 and it provides data for 55 of our 71 fields. We use the 
WISE 3.4//m, and 4.6//m data to look for signs of enhanced 
mid-infrared scattering and the WISE 12.0/zm and 22.0;um data 
to characterise the mid-infrared dust emission and to look for 
indications of ongoing star formation. The data were converted 
to surface brightness units w ith the conv ersion factors given in 
the explanatory supplement (ICutri et al.l l201 1). The calibration 
uncertainty is ~6% for the 22 yum band and less for the shorter 
wavelengths. 



We calculated dust extinction maps for each field wit h stars 
from the Two Micron All Sky S urvey (2MASS , Skrutsk ie et alj 
12006 ). using the NICER method dLombardi & Alves,i200ll) . The 
values of Ay were obtained from near-infrared (NIR) colour ex- 
cess measurements assuming an extinction law with Ry=3.l. 
With the assumption of Rv=5.5, possibly more appropriate for 
the densest regions, the Ay estimates would increase by ~16% 
(see lJuvela et al.ir20I Ih . The spatial resolution of all the pro- 
duced extinction maps is two arc minutes. For distant sources, 
the extinction of the target clouds cannot be reliably repro- 
duced because of the poor resolution and the increasing num- 
ber of foreground stars. The same applies to some extent to the 
nearby, high latitude targets. No special steps have been taken 
to eliminate the contamination by foreground stars (see, e.g., 
ISchneider et al.ll20I ll) . apart from the sigma clipping procedure 
(clipping performed at the 3-cr level) included in the NICER 
method. Therefore, the extinction maps should not be trusted 
unreservedly as tracers of the column density of the examined 
clouds. On the contrary, the possible discrepancy between the 
Ay map and the sub-millimetre emission serves as an indicator 
of a large source distance. 



3. Results 

3.1. Cloud properties derived from SPIRE 

In this section we present the results derived from the SPIRE 
data. These include colour temperature and column density maps 
and, when distance estimates are available, estimates of the 
cloud masses. In Sect. 13.1.41 we examine the properties of the 
filamentary cloud structures. 



3.1.1. Colour temperature maps 

The colour temperature maps of the large grain emission are esti- 
mated using the 250 fj.m, 350 jum, and 500 fim SPIRE maps with 
the main goal of identifying the coldest clumps. The data cover 
only the longer wavelength side of the dust emission peak which 
means that the derived temperatures may not be very accurate, 
especially for warm regions with temperatures close to 20 K or 
above. On the other hand, the absence of shorter wavelengths 
decreases the bias that results from temperature variations along 
the line of sight ( Shetty et al. 2009b a; Malinen et al. 201 1 ). The 
statistical noise will be larger than, for example, if one included 
the PACS bands, but the values will be more representative of 
the bulk of the cold dust. The absence of data below or near 
lOOjum also avoids the problem of a possible contribution from 
stochastically heated small grains. 

The maps were convolved to a 40 arcsec resolution and, 
for each pixel, the SED was fitted with ByiTduat)'^ keeping 
the spectral index y6 at a fixed value of 2.0. Several studies 
have suggested that the spect ral index may increase in cold 
and dens e environments (e.g. Dupac et aP l2003t iDesert et 
2008; Pl anck Collaboration et alj20I Id) . If this anticorrelation 
between the dust temperature and the spectral index is true 
also in our fields, the derived temperature maps will under- 
estimate the range of temperature variations and, in particu- 
lar, will overestimate the temperature of the coldest regions. 
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Fig. 3. Uncertainty of the colour temperatures estimated with a 
Monte Carlo study. The curves show, for each field, the mean er- 
ror of the temperature if the uncertainty of the surface brightness 
data were 13%. The horizontal line is the 1-cr value (calculated 
over all fields) of the bias associated with the uncertainty of the 
intensity scale zero points. 



environments (iBeckwith et al.l[r990l) . The value of dust opac- 
ity is uncertain and is observed to vary from region to re- 
gion (e^. Kramer et al. 2003; Stepnik et al. 2003; Lehtin en et al.l 
|2004 120071: iRidderstad et al J 120061: iMartin et al.li201lb . In dif- 
fuse medium the expected value of k is lower by a factor of 
two jBoulanger et al. Il996l) . However, the value chosen here is 
close to the predictions and the observations of dense clouds 
(e.g. Ossenkopf & Hennii ig 1994 ; Nu tter et_al. 2006 2008) a nd 
also is the value used in IPlanck Collaboration et al. 
The column density maps of the first eleven sources are shown 
in Fig. |6]and for the rest in Appendix iBl (online edition). The 
resolution of the maps is 40". 

Mass estimates are possible for the sources with distance es- 
timates, i.e., 59 out of the 71 fields (see Table [1]). Masses are 
calculated for the entire field, for some major filaments (see 
Sect. I3.1.4I I. and for one or two representative clumps within 
each field (see Sect. 13.2.11 ). The mass estimates are listed in 
Table|2l These mass estimates are intended only for general char- 
acterization of the target fields, while more detailed analysis of 
all clumps, including their mass spectra, will be presented in a 
later publication. 



leading to an underestimation of the masses of these regions 
(Planck Collaboration et al. 2011c). However, the temperature 
maps will still fulfill their main purpose, identifying the major 
relative temperature variations within the regions. 

The effect of statistical errors was examined with the help 
of Monte-Carlo simulations where a 13% uncertainty in the sur- 
face brightness values was assumed. Based on the correlations 
between the diff'erent bands, this is a conservative estimate and 
probably even twice the value of the typical true noise. The 13% 
uncertainty of the surface brightness values would translate to 
an error below 1 K in cold regions, the error increasing to ~3 K 
at 20 K (Fig. [3]l. The calibration accuracy of the SPIRE data is 
estimated to be 7%. This is well within the above limits. 

The effect of the intensity zero points was estimated in a 
similar fashion, examining realisations of the temperature maps 
when the zero points were modified in accordance with the 
uncertainty of the comparison of Herschel and Planck/IRIS 
data. This uncertainty was less significant, typically less than 
0.5 K and ~1K in a few cases (G130.42-47.07, G161.55-9.30, 
G37.49+3.03). However, the effects on the relative temperature 
variations within a given field are, of course, much smaller The 
temperatures can be further affected by mapping artifacts, e.g., 
residual striping. The striping is usually visible only towards the 
map edges where its effect can rise to ~0.5 K. 

Two examples of calculated temperature maps are seen in 
frame a of Figs. |4] and |5] The colour temperature maps of all the 
fields are presented in Appendix|A] (online edition). 

3.1.2. Column density maps and masses 

The column densities averaged over a 40" beam are calculated 
using the formula 



A?(H2) 



By{T)K^imn ' 



(1) 



with the intensity and temperature of the previous SED fits 
(Sect, [mi l and using ju = 2.33 for the particle mass per hy- 
drogen molecule. We use a dust opacity k obtained from the for- 
mula 0.1 cm^/g (v/lOOOGHz)'^ that is applicable to high density 



3.1 .3. Categorisation of fine source fields 

The sources, as seen through the sub-millimetre emission, ex- 
hibit different morphologies from isolated unstructured clumps 
to more complex configurations. We grouped the sources to non 
exclusive sets of cometary, filamentary, isolated, and complex 
sources (tags C, F, I, and X in Table[T]l. 

The category of cometary sources refers to the shape of the 
major clump or clumps that suggests deformation by an external 
press ure that has resulted in a shape typical of cometary glob- 
ules dHawarden & Brandll 19761) (e.g. LDN 1780 in Fig. ??, on- 
line edition). We distinguish sources that show more local indi- 
cation of direct dynamic interaction, e.g., in the form of a sharp 
cloud boundary, without the overall cloud shape being cometary 
(tag B\ e.g. field G25.86-h6.22 shown in Fig. ??, online edition). 
About one field in ten shows clear signs of either type of a dy- 
namical interaction. 

In the group marked as filamentary (tag F; e.g. field G82.65- 
2.00 in Fig. ?? in online edition), the clump or clumps are lo- 
cated inside narrow, elongated structures or the filaments them- 
selves are the dominant source of dust emission. The filaments 
will be discussed in more detail below but one must already note 
that the definition of a filamentary cloud is not clear and, as a 
result of turbulent motions, some degree of filamentary structure 
is visible in all fields. In about half of the fields the overall mor- 
phology is dominated by a single or a couple of filaments. 

We denote fields as complex when they contain several 
clumps that are not clearly arranged onto a small number of fila- 
ments (tag X, e.g. field G86.97-4.06 in Fig. ?? in online edition). 
At low latitudes, this also could be caused by line-of-sight confu- 
sion without it being an intrinsic property of the clouds observed. 
On the other hand, the isolated sources (tag /) are typically found 
at intermediate and high latitudes. 

Few fields can be assigned unambiguously to a single cate- 
gory and in most cases one can identify characteristics of more 
than one category as indicated in the eighth column of Table [T] 
The classification is, of course, a subjective one and is made only 
to help the discussion below. 
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Fig. 4. Selected observed and derived quantities for the field G4. 18+35.79 (LDN 134). Shown are the colour temperature map 
derived from SPIRE data (frame a, 40" resolution), the 250 /im SPIRE surface brightness map (frame b, 18" resolution), the 
AKARI wide filter maps at 140 |/m and 90 /im (frames c and e), the visual extinction Ay derived from 2MASS catalog stars (frame 
d), and the WISE 22 /vm intensity (frame /). For the other fields, see Appendix lAl The respective beam sizes are indicated in the 
lower right hand corner of each frame. The positions of selected clumps (black circles; see Sect. 13.2. Il l and the reference regions 
used for background subtraction (white circle; see Sect. 13.2.1! ) are also shown. In frame b, the arrow indicates the position of the 
stripe plotted in Fig.|9] 



3.1.4. Properties of the filaments 

Half of the examined fields are categorised as having a promi- 
nently filamentary structure (see Sect. l3.1.3] and Table[T]i. We ex- 
amine here the properties of selected cloud segments. The selec- 
tion is carried out using column density maps that are obtained 
by combining, without further convolution, the 250 yum surface 
brightness data (resolution ~ 18") and the colour temperature 
estimates (resolution 40", see Sect. 13.1. 11 1. This gives a nomi- 
nal resolution of ~ 20" although, of course, the column density 
estimates depend on the temperature information that is avail- 
able only at half of this resolution. As the dust temperature is 
expected to decrease towards the centre of the filaments, the low 
resolution of the temperature data is likely to cause the central 
column densities to be underestimated and the column density 
profiles to appear flatter than in reality. However, the colour tem- 
peratures are known to be b iased estimates of th e mass weighted 
average dust temperature ( Shettv et al.l l2009b a; M alinen et al] 
1201 1 ) and the low resolution of the temperature data may not be 
the main source of uncertainty. We will return to this question in 
Sect.m 

We select from each of the maps the most prominent fila- 
ment (or, more generally, elongated structure) and trace by eye 
the path along its central ridge. The column density profiles per- 



pendicular to the path are extracted at 20" intervals. We examine 
the average column density profile as well as its variations along 
the length of the filament. The estimated width of the filaments 
depends on how far from the centre of the filament the profile 
is followed and how the background is subtracted. If the area is 
fixed in angular units, one will probe larger linear scales for the 
more distant sources and, because of the hierarchical structure 
of the ISM, the result will be a correlation between the distance 
and the width estimate. To avoid this, we use a fixed range of 
[-0.4pc,+0.4pc] around the peak of the column density profile. 
This scale is adequate for most fields as it is resolved for the most 
distant sources while remaining smaller than the map size in the 
case of the more nearby targets. The FWHM is measured after 
the subtraction of a constant background level that is estimated 
as the minimum value within the [-0.4 pc, -1-0.4 pc] interval. 

We also fit the filament profiles with Plummer-like density 
profiles 



Pp(r) = 



Pc 



N(r)^Ap- 



(2) 



[|1 + (r/Rn.d^\P'^] [1 + (r/iJflJ^lC-i)/^ 

( see iNutter et al.1 120081; lArzoumanian et al.ll201 ll) . In the equa- 
tion Pc is the central density, /?flat is the size of the inner flat 
part, and p describes the steepness of the profile. For an isother- 
mal cylinder in hydrostatic equilibrium, the value of p should be 
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Fig. 5. Selected observed and derived quantities for the field Gl 30. 37+ 11. 26 (see Fig.|4]for details). 



4 (IOstrikeiill964l) but the values derived from observations are 
usually smaller (e.g. lArzoumanian et al The factor A,, is 

obtained from the formula A „ = f '^"-,„n , which depends 

P cost J-OO (1+H-)W-' r 

on the unknown inclinati on of the filament for wh ich we will 
assume a value of ; = 0. lArzoumanian et a noted that 

while this assumption does not affect the analysis of the shapes 
of the profiles, the observed column density is on the average 
~60% higher than the column density measured perpendicular 
to the filament. This could affect the estimates of the stability of 
the filaments. On the other hand, it is clear that selection effects 
here favour small inclination angles. The data are fitted with the 
Plummer profile together with a linear background, keeping the 
centre position along the filament profile as a free parameter. 
Because a linear background is part of the fit, these results do 
not depend on the absolute level of the background but are still 
affected by the extent of the fitted region, [-0.4 pc, +0.4 pc]. 

When we calculate the average profile of the whole filament 
segment, the individual profiles are first aligned so that the peak 
column density appears at the centre of the profile. This decision 
tends to minimise the width of the average profile. 

Fig.|7]shows an example of the results for the field G163.82- 
8.44. The frame a shows, on the column density map, the initial 
path that was selected by eye (the white line). The beginning of 
the path is marked with a filled circle and the end with an open 
circle. The yellow line traces the peaks of the individual column 
density profiles. 

Figure [TJ? shows the mean column density profile and the 
fitted Plummer profile. The profile was obtained by averaging 
the individual profiles where the column density rises above the 
median value over all the profiles that were sampled at 20" in- 



tervals. We also ignore those profiles where the column density 
peak has shifted more than 5' from the initial filament path. In 
the same frame are drawn the individual profiles where the peak 
is amongst the 20% of the highest column densities along the 
filament (the yellow lines). 

The frame c shows the filament FWHM and the parame- 
ter 7?flat of the Plummer profile fit as a function of the filament 
length. For comparison, the column density along the filament 
is shown as a solid line. Both the frames b and c indicate varia- 
tions in the width of the filament. These often are anticorrelated 
with the column density (and, consequently, temperature) and 
are discussed further in Sect. 14.31 

Figure|8]shows the corresponding results for the more nearby 
field of G300.86-9.00 (PCC550). For the other fields the figures 
are included in the Appendix ICl 

Table [3] lists the properties derived for the average filament 
profile of each field where distance estimates were available. 
These include the average column density (corresponding to the 
average profile in Fig.lTJi), the FWHM values determined from 
data within 0.4 pc of the filament centre, and the parameters TJflat 
and p from the fits of the Plummer profile. The figures show that 
there is often considerable variation in the width along the length 
of the selected structures. Therefore, the listed widths are not 
always representative of the cleanest filament segments which 
tend to be more narrow. The table includes further estimates of 
the mass per linear distance that can be compared to the crit- 
ical value of a gravitationally stable filament, - 2cl/G 
(Ostriker 1964). Instead of using the estimated dust temperature, 
we use a sound speed for a constant gas temperature of 12K 
that gives Mj^"' = 20 Mq pc" ' . This is justified because, for most 
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Fig. 6. Column density maps of the 12 first fields in the order of increasing Galactic longitude. The map resolution is 40" and the 
areas are the same as in Fig. |4] and Appendix lAl (online edition). The colour scale has a square root scale. The column densities of 
the other fields are shown in Appendix IB] 



of the filament volume, the gas and dust temperatures remain un- 
coupled ( Goldsmith 2001; Juvela & Ysard 2011 ). Note that this 
Mj™^ is applicable only in an isothermal case, in the absence of 
other supporting forces like magnetic fields and turbulence. Of 
the 30 filaments listed in the table, 17 appear to be supercritical. 



3.2. Comparison with otiner infrared data 

We have compared the sub-millimetre observations with the 
mid-infrared and far-infrared data from the WISE and AKARI 
satellites. The shorter wavelengths are more sensitive to tem- 
perature differences in the thin surface layers of the clouds and 
can trace variations in the heating radiation field or in the rela- 
tive abundance of small and large grains (Laureiis et al.lll989l: 
iBernard et all Il992t iRidderstad & Juvelal .2010) . Therefore, 
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Fig. 7. Properties of a filament in the field G 163. 82-8. 44. The frame a shows the column density map. The white line shows the 
filament that was originally traced by eye and the yellow line follows the ridge that is formed by the peaks of the column density 
profiles in the perpendicular direction. The black filled circle indicates the start of the examined filament section. The tick marks 
are drawn at 5 arcmin intervals. The frame b shows the average column density profile of the filament (black line), as well as the 
Plummer profile (the red dashed line on top of the black line) that was fitted together with a linear baseline (the blue dashed line) 
over the range -0.4 pc to -1-0.4 pc. The yellow lines show individual column density profiles for 20% of the cuts with the highest 
column densities (values normalized to a peak value of one, the right hand scale). The frame c shows the FWHM values (the black 
circles) and the parameter 7?flat of the Plummer fit (the red crosses), and the column density along the ridge of the filament (the solid 
line and the right hand scale) as a function of the distance along the filament. With the estimated distance of 350 pc, the SPIRE 
250^(m data resolution of 18" corresponds to 0.03 pc. The plots for the other fields are shown in AppendixICl 




Fig. 8. The filament in the field G300.86-9.00 (PCC550) that is part of the Musca molecular cloud at the distance of 230 pc. See 
Fig.|2]for the details. 
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these data can give further clues to the nature and the formation 
process of the clumps. 

In addition to the SPIRE 250jum data. Fig. Hand the fig- 
ures in Appendix |A] show the WISE 22 fim and the AKARJ 
90jum, and 140 yum surface brightness maps. For example, the 
field G4. 18+35.79 is apparently quiescent with no mid-infrared 
(MIR) or far-infrared (FIR) point sources (Fig.|4|. On the other 
hand, in 0130. 37+1 1.26 (Fig.|5]) the 22fj.m data indicate clear 
star formation activity in connection with this mainly cold 
clump. The enhanced 22fim emission towards the south further 
suggests an anisotropy of the radiation field. The lack of corre- 
lation or even clear anticorrelation between the PAH emission 
dominated MIR and the FIR/sub-mm emission is a feature re- 
peated in many of the fields. 

Figure |9] and Figs. ??-?? (online edition) show cross sec- 
tions of the surface brightness data along the lines marked in 
frame b of Fig. |4]and the figures of Appendix lAl The figures 
show the intensities after the subtraction of the values in the ref- 
erence regions that were marked in the Fig. |4] and the figures of 
AppendixlAl 

3.2.1 . Clump properties from dust emission 

In this section we look in more detail at the properties of some 
major clumps, combining Herschel observations with the in- 
frared data. The analysis uses the apertures that were marked in 
the frames a of Fig.|4]and figures in AppendixlAl(online edition). 

The masses within the apertures were first determined using 
SPIRE data only. We derive two mass estimates. The first one is 
obtained by directly integrating over the aperture in the column 
density maps of Sect. l3.1.2l The spatial resolution of the column 
density maps is 40". Because the column densities are calcu- 
lated using the absolute surface brightness values, these should 
correspond to the total mass along the line of sight. However, the 
estimates are based on colour temperatures that, in the case of 
cold cores, can seriously overestimate the temperature inside the 
cores. We calculate another mass estimate using aperture pho- 
tometry. The SPIRE surface brightness maps are convolved to 
a resolution of 40". We measure the signal within the aperture 
and subtract the background using an annulus that extends to 
1.5-2.0 times the radius of the aperture, thus removing most 
of the diffuse signal. The aperture sizes are listed in Table |2l 
The uncertainties of the flux values are derived from the sur- 
face brightness fluctuation in the annulus. The resulting SPIRE 
fluxes are fitted with a modified blackbody to estimate the clump 
masses. Because of the local background subtraction, these esti- 
mates should better reflect the temperature and the mass of the 
dense structures only. The mass estimates are listed in Table |2] 

We have constructed the clump SEDs from 22/im to the 
SPIRE wavelengths. Figure [TO] shows the SEDs for clumps in 
the first twelve fields (in order of increasing Oalactic longi- 
tude) and similar figures for the remaining fields can be found in 
Appendix iP] In addition to the three SPIRE channels (250 /^m, 
350/im, and 500 yum) the plots include the AKARI FIS bands 
(65 /im, 90 /im, 140 /zm, and 160 jum) and the WISE data at 12 
and 22/im. The data were convolved to one arcmin resolution. 
The local background subtraction was performed as in the case 
of SPIRE data above, using an annulus around the aperture. 
Alternatively, because we do not have absolute surface bright- 
ness measurements in all infrared bands, the background can be 
estimated using the reference areas indicated in the frames b of 
Fig- Sand in Appendix |A] If the signal in the reference area is 
small, the resulting spectra should give a good estimate of the in- 
tensity and the spectral shape of the total emission along the line 
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Fig. 11. Comparison of the temperature estimates for the se- 
lected clumps. In the left hand frames, the x-axis is the tempera- 
ture derived from the total surface brightness. This is correlated 
with the minimum values within the apertures (frame b) and the 
values obtained with a local background subtraction (frame a). 
In the right hand frames the SPIRE temperatures (with back- 
ground subtraction) are compared with the values obtained using 
the combination of SPIRE and AKARI data (A > lOOyum) and 
with (frame c) or without (frame d) the subtraction of the local 
background. 



of sight. The data at wavelengths A > 100//m were fitted with 
modified blackbodies, By{T) v^, with a fixed value of the emis- 
sivity spectral index, yS = 2.0. The resulting colour temperatures 
are given in the figures. 

We now have several colour temperature estimates for each 
aperture. Two are based on SPIRE only (as in Fig.|4]and Table|2]) 
and two additionally include the AKARI 140/im and 160yum 
data (in Fig. [10] and Appendix iDl). Furthermore, in two cases we 
use the total surface brightness values (or an approximation ob- 
tained with the help of the reference areas) while in the other 
two cases the local background is subtracted using an annulus. 
The temperature estimates are compared in Fig.[TT] The subtrac- 
tion of the local background typically leads to ~2 K lower colour 
temperatures as expected for sources colder than their environ- 
ment. The values are lowered by a further 1-2 K when the FIR 
data are included in the modified blackbody fits. One would have 
expected the temperatures to be lowe r when the shorter wave- 
lengths are excluded (see IShettv et al. 2009b). The result could 
suggest a small difference in the calibration of the two instru- 
ments. The other possibility is that the FIR intensities are partly 
saturated. For the typical beam averaged peak column densities 
of 10^^ cm"^, the 100 jum optical depth is of the order of 0.01, 
too little to cause noticeable effects. If the mass is concentrated 
to a very small area within the beam, the optical depths will be 
higher A 1-2 K difference in the colour temperature would re- 
quire column densities in excess of 10^"* cm"-. Such values are 
possible, for example, in case of distant, pre-high-mass clumps. 
However, the optical depths are unlikely to be the main explana- 
tion for the whole sample because that would require that most 
of the mass within the beam is always concentrated on such high 
optical depth sightlines. 



3.2.2. Searcli for coresliine 



ISteinacker et al ] (l2010h detected enhanced mid-infrared surface 
brightness in the Spitzer/IRAC 3.6 and 4.5 fim channels towards 
the cloud LDN 183. The high intensity of the 3.6 fim band rel- 
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Fig. 9. Cross sections of the surface brightness data along the lines indicated in Figs.|4]-(5]and in the figures of AppendixlAl The lines 
show the 250 yum SPIRE data (thick black line), AKARI 140 fj.m and 90 /im data (solid and dashed blue lines), and, when available, 
the WISE 22 /im data (dotted line). The average surface brightness in the reference region (see Fig.|^) has been subtracted from 
the plotted values. The red thick line is the colour temperature. The data have been convolved to the resolution of one arc minute. 
The 90 yum data have been scaled by a factor 20 and the 22 yum data by a factor of 600. When a different scaling has been used, 
the wavelength and the multiplicative scaling factor are given in the frame (/l:factor). The plots for the other fields are shown in 
AppendixiBl 



ative to the 4.5 fim band, the absence of emission in the 5.8 yum 
band, the presence of absorption in the 8 yum band, and other 
considerations, led them to the conclusion that neither PAHs nor 
warm grains could explain these observations which were in- 
terpreted as a sign of enhanced light scattering caused by an 
increase in the size of the dust grains. This excess, named by 
the authors as coreshine, would thus serve as a tracer of the 
grain evolution, similar to the increase of the sub-millimetre 



dust opacity (e.g.'Stepnik et al.l2003LlPianck Collaboration et al.l 
2011a). Thus, this phenomenon could also be used as an in- 
dicator of the evolutionary stage of our clumps. We examined 
the correlation between the WISE satellite 3.4, 4.6, and 12yum 
bands and the 250yum band measurements from SPIRE. If the 
grain size has indeed significantly increased, the WISE 3.4 yum 
band should include an additional component of scattered light 
that makes it brighter than expected compared to its surround- 
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Fig. 10. Spectral energy distributions for the apertures indicated in Fig.|4]and the Appendix lAl The plots include the 22 /im WISE 
data, the AKARI data at 65yum, 90yum, 140 //m, and 160 yum, and the three SPIRE channels at 250yum, 350yum, and 500 yum. In 
most fields two aperture positions were chosen and the data for the second one are shown in red. The background subtraction was 
done using the reference areas marked in Figs. |4] and the figures of Appendix |A](solid symbols) or by using a local annulus (open 
symbols). The colour temperatures from the modified blackbody fits with /3 - 2 are listed in the frames. Plots for the other fields 
can be found in AppendixiDl 



ings and this should be correlated with the absence of emission, 
or better the absorption in the WISE 12 /urn band, weaker or no 
emission in the WISE 4.6 fim band and these features should ap- 
pear inside the 250 /urn emission region. 

T he 4.6 ^m band includes H 2 lines that can be strong in out- 
flows dCvganowski et al.l201 ih . Also the 3.4//mbandis affected 
by H2 emission but, for typical excitation conditions, the inten- 
sity is much lower than in the 4.6 fim band jReach et al. 2006). 
Furthermore, the H2 emission is also seen in the next MIR bands 
(5.8 jum and 8 fim for Spitzer/IRAC and 12 /urn for WISE) and is 
therefore easy to distinguish from coreshine for which absorp- 
tion is expected at these longer wavelengths. 



One caveat concerns the absorption of the MIR radiation 
from the sky behind the source. Because the cloud opacity in- 
creases towards the shorter wavelengths, the background radi- 
ation is attenuated more at 3Afj.m than at 4.6 yum. If the back- 
ground level is high, the detection of the coreshine may require 
modelling of the absorption and emission processes. This sug- 
gests that the direct effect is best seen at high Galactic latitudes 
where the sky brig htness is low. This has been remarked by 
Pag ani et alJ (|2010^ who could find no cases of coreshine emis- 
sion in the galactic plane itself. An upper limit of 0.37 MJy sr ' 
in the background brightness level (which is linked to the den- 
sity of the field stars) to allow for coreshine detection has 
been derived and will be discussed elsewhere (Pagani et al. in 
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preparation). Another limitation comes from the low resolution 
of the WISE satellite (6.1" for all three bands, 3.4, 4.6 and 
12yL(m), which, combined with a high density of stars close to 
the galactic plane, tends to fill up the space with starlight, hid- 
ing possible coreshine effects. Figure [TSlshows data for the field 
G6. 03+36. 73 (LDN183), one of the strongest coreshine detec- 
tions, already discussed by Steinacker et al. (2010). The central 
filament, as traced by the 250 fim emission, is clearly visible in 
the 3.4/im map. There is clear detection in emission at 4.6/im 
and in absorption at I2ij.m . The coreshine intensity is compara- 
ble to the one reported by ISteinacker et ail (l2010h . The 3.4yum 
signal increase is larger than for the 4.6;um one and the cross- 
sections of the surface brightness show that the increase of the 
3.4 jum/4.6 fim ratio is correlated with both the 250 fim emission 
and the decrease of the dust temperature. 

We have similarly examined the 3.4, 4.6 and lljjm images 
of all the 56 fields for which public WISE data are available. We 
compare them with the 250 yum SPIRE surface brightness maps 
that act as a tracer of the large grains. Figure [13] shows another 
example of detected coreshine. 



4. Discussion 

4.1. The nature of the observed fields 

The observations show that the clouds associated with cold 
clumps have varied morphologies. The cloud sizes range from 
nearby globules ~0.2pc to massive filaments more than 20 pc 
in length. With very few exceptions, the clumps exhibit a sig- 
nificant amount of sub-structure down to the resolution limit 
(0.01 pc for the nearest regions, over 0.2 pc for the most distant 
ones). 

Of the 71 fields examined, 36 fields or 51% were catego- 
rized as having a clear filamentary structure. The definition of 
a 'filament' here applies to the major structures but is taken 
in a widest possible sense, meaning anything from slightly 
elongated structures (e.g., Gl. 94+6.07 or GlO.20+2.39) to the 
very long and narrow 'real' filaments like those observed in 
the fields G82.65-2.00 and G276.78+1.75. Several filaments 
appear to have formed, or at least significantly deformed, by 
some external force. This is visible as an asymmetry or a sharp 
drop in the spatial distribution of emission and column density 
(e.g., G25. 86+6.22, G94. 15+6.50, G130.37+1 1.26, or G315.88- 
21.44, see Figs. ??, ??,|5j and ?? in online edition). Thus, the ori- 
gin of these filamentary structures could be different, i.e., turbu- 
lence vs. direct compression by radiation or pressure waves. This 
could be reflected in their properties which are further discussed 
in Sect. 14.31 However, further study of the interaction with, e.g., 
nearby HII regions is deferred to a future paper. 

Altogether eight fields showed features reminiscent of com- 
pression by an external force (marked with 'B' in column 8 of 
Table [T]l that has lead to a sharp boundary visible in the col- 
umn density maps. Such an interaction can also manifest itself 
at larger scales, as a cometary of a filamentary shape of the entire 
cloud. This sample includes some well known cometary shaped 
clouds (e.g., LDN 1780). At low resolution, the cloud LDN1340 
(field GOO. 37+1 1.26) also looks like a classic cometary cloud 
but at higher resolution is found to consist of an intricate network 
of filaments that, particularly in the tail of the cloud, breaks down 
to many smaller clumps. 

The clumps with a simple structure are mostly nearby ob- 
jects where the spatial resolution allows us to resolve the core 
scales. However, there also are some more distant objects in this 
category. The target G39.65+1.75 is at a distance of 1.8kpc and 



it is probable that it would not be such a featureless 'blob' if 
it could be better resolved. On the other hand, the main clump 
in G26.34+8.65 is resolved but shows remarkably little struc- 
ture for such a distant and therefore large object. However, one 
must keep in mind that the distance estimates of some of the 
objects are still quite uncertain and the distances of objects like 
G26.34+8.65 might be overestimated. 

The colour temperature maps (Figs. I^ll] and Appendix lAb 
show that the dense structures are, in colour temperature, typi- 
cally 2-3 degrees colder than their diffuse environment. The ap- 
pearance of the temperature maps, i.e. their strongly spatially 
correlated values, indicates that the maps are dominated by real 
temperature variations rather than noise. The analysis of the un- 
certainties suggested that in the temperature maps the relative 
values are quite reliable (Sect. [TO]) although calibration errors 
could cause a shift of about +1 K (see Sect. l3.2.11 i. 

The colour temperature calculations were based on the 
use of a fixed value of the dust emissivity spectral index, 
/3. There is evidence that, at least statistically, there ex- 
ists an inverse relation betw een the ob served spect ral in- 
dex and the dust t emperature dDupacet al. 2003; Desert et alJ 
2008; Planck Cofl aboration et all 1201 Id IVeneziani et al.1 120101: 
Paradis et al., ,2010i) . It is still debated to what extent this is 



caused by the intrinsi c dust properties, by the line-of-sight 
temperature variations (Shettv et al. 2009b: iMahnen et"al]|201 U 
Juvela & Ysard 2012a; Ysard et al. 2012), and, in particular, by 
the noise that is known to produce some anticor relation be- 
tween the colour tem perature and the spectral index ( Shettv et al.l 
l2009aL ITuvela & Y sard 2012b). However, the anticorrelation ap- 
pears to be an established observational fact and s tronger than 
what i s produced by the noise only (P lanck Collaboration et al.l 
1201 Id) . If the average spectral index is larger in cold regions, the 
temperature contrasts (i.e., the difference in derived colour tem- 
perature) would be higher than indicated by our maps. Figure[T4l 
shows some spectra consistent with the relation /3{T) = (6 + 
u)Ty^ with parameters S - 0.02 and o) - 0.035 as esti- 
mated from the analysis of Planck observations of dense cores 
(iPlanck Collaboration et al.l l201 id) . The relative intensities in 
the SPIRE channels and the SED fits obtained with a fixed value 
of P - 2.0 are shown. For example, when a value of p - 2.0 is 
assumed, the estimated colour temperature is 12.2 K while the 
actual values of the spectral index and temperature are fi-l.l 
and Tdust^lOK. The figure also illustrates that, without further 
far-infrared data, simultaneous determination of colour temper- 
ature and spectral index is not possible. 

We conclude that the colour temperature maps do reliably 
identify the coldest clumps although the absolute temperature 
values may be uncertain up to a couple of degrees. The real 
dust temperature inside the clumps cannot be directly measured 
and could be estimated only by modelling. From earlier theo- 
retical studies, it is clear that in many cases the cent ral tem- 
perature of th e de nse co res goes below 10 K (e.g., Eva ns et aH 
2001'; 'Harju et al." 2008J). In partic ular, in LDN 183, our field 
G6. 03+36. 73, Pagani et al.l (|2004|) reported dust temperatures 
down to ~7K. In the dense and cold environment the dust 
grains are expected to undergo coagulation that can lead to in- 
c reased opac ity at sub- millimet re wavelengths (Ossenkopf 19931; 
iStennik et al.ll2003UYsard et al.ll2012tlKohler & et al.ll2012l) . At 
the same time, the abundance of small grains should decrease, 
leading to lower mid-infrared emission. Examination of the 
data for the apertures shown in the figures of Appendix [A] re- 
vealed no significant anticorrelation between the surface bright- 
ness ratio 22yum/90yL(m and the ratio 90//m/250yL(m or the 
colour temperature measured in the wavelength range 250- 
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Fig. 12. Evidence of coreshine in the field G6.03+36.73 (LDN 183). Shown are the WISE 3.4 fim, 4.6 fim, and 12 /urn WISE maps 
(frames a-c) and the SPIRE 250/im map (frame d). Frame e shows surface brightness profiles along the arrow marked in the 
previous frames. The lines correspond to the 3.4 fim (solid line), 4.6 /urn (dashed line), 12 /urn (dash-dotted line), and 250/zm data 
(dotted line). The numbers in the upper right corner give the possible scaling applied to the data before plotting (A:factor). The red 
line and the right hand scale show the colour temperature. 




500 //m. There could be several reasons for this. Firstly, most 
clumps could be too young so that the coagulation process, 
which may take millions of year s, has not had time to alter 
the abundance of the small g rains (lOssenkopf & Hennindll994t 
IChakrabarti & McKedllK)! . Another natural reason is that the 
22 jum emission comes mainly from a region with Ay < 2™, as 
indicated by MIR li mb brighte ning in radiative transfer models 
dBernard et al.lll992l: (pischera & Dopital l2008). and thus is not 
sensitive to processes deep inside the cores. The IR point sources 
associated with the cold clumps are another complication. Many 
regions are already forming stars so that the MIR data, the 22 yum 
band included, are dominated by the warm sources rather than 



the extended dust emission. This contaminates the correlations 
between MIR emission and the mean colour temperature of the 
large grains. 



In a few individual clumps, an anticorrelation does exist be- 
tween the column density, as traced by the sub-millimetre emis- 
sion, and the MIR emission. These include the fields G6.0+36.7, 
G37.49+3.03, and G26.34-h8.65 as shown in Fig. |9] and some 
further examples can be seen in AppendixlBl (online edition). 
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Fig. 14. Spectral energy distributions (solid curves) for 10, 13, 
16, and 19 K colo ur temperatures assuming the /3{T) relation 
given in iPlanckC ollaboration et al. (201 Ic). The corresponding 
values at the SPIRE wavelengths are plotted with 10% error bars. 
The dashed lines show the spectra fitted to these data points 
assuming /3 - 2.0. The corresponding temperature values are 
quoted after the arrows. 



4.2. Connection witli star formation 

Star formation activity is visible in the MIR maps as sources 
spatially correlated to the cold dust emission. The best example 
is the field G216. 76-2.58 where ten strong WISE MIR sources 
are detected along the high column density structures. This does 
not prevent cold dust being detected at longer wavelengths in 
exactly the same areas, with colour temperatures down to Tdust ~ 
12 K. As a preliminary quantitative estimate of the star formation 
activity we counted the surface density of MIR sources. For each 
field for which WISE data were available, we made a list of the 
22 fim point sources. The detection was done with the SExtractor 
program, excluding sources that are significantly more extended 
than the beam (FWHM larger than 10 pixels - 60"). Each field 
was divided into a low and a high column density part using 
the median column density as the dividing value. The source 
densities were calculated for the two parts. Figure [15] shows the 
source density in the high column density areas (i.e., sources 
likely to be associated with the clumps) relative to the average 
source density in the field. 

The interpretation of such correlations is not straightforward. 
The youngest objects should reside within the areas of high col- 
umn density and thus lead to a positive correlation (fields on the 
right hand side of Fig. \T5[ . A negative correlation could appear 
in more evolved regions where the sources have cleared their sur- 
roundings. However, with the typical clump column densities of 
10^^ cm"^ (Ay ~ 10), the 22yum optical depth approaches unity 
and this is enough to affect the number density of the detected 
background sources. Therefore, the fields in the left hand part 
of Fig. [15] could represent either cases without associated star 
formation or cases with a larger number of background sources. 
The number of MIR sources and the correlations may be further 
affected by the source distances because of the sensitivity limit 
of the WISE survey. However, the overall source densities of the 
fields and the fraction of sources in the high column density ar- 
eas do not appear to be correlated (see Fig.[T5]l. 

About 50 percen t of the clumps are located on FIR loops 
iKonvves et alj (l2007l) that are boundary regions (shells) around 
the voids in the galactic ISM. A preliminary analysis of a sam- 



ple of sources and their AKARI FIR photometry showed that 
the clumps located on the shells have ~20% higher probabil- 
ity of being associated with YSO candidates (Marton 201^. 
In particular, a detailed analysis of the 62 IR point sources 
in the "snake" cloud (California nebula, G163.82-8.44) region 
resulted in a number of YSO candidates . The sources were 
analysed using the AKARI IRC and FIS (Ishi hara et all 120101: 
Murakami et al. 2007; Yamamura et al. 2010), the 2MASS and 
IRAS point source catalogues, as well as optical photometry data 
available in public archives. On the basis of the physical parame- 
ters derived for the sources, there are 1 1 low-mass and 2 interme- 
diate mass young stellar objects in evolutionary phases ranging 
from Class O/I to Class III. In that particular field, the estimated 
ages suggest that the star formation started several million years 
ago and is still on-going. Further details will be given in a sepa- 
rate paper (Zahorecz et al. 201 1 in preparation). 

4.3. Filamentary structures 

4.3.1 . The uncertainties of filament analysis 

The characterisation of the filaments depends on the way the 
analysed structures are selected and on the analysis methods 
themselves. We selected at most one elongated structure from 
each field. One should consult the figures in Appendix |C] to 
check to what extent these correspond to ones own conception of 
filaments. There are automatic routines for the detection o f fila- 
m ents, for example the DisP erSE routine (|Sousbidl201 ih used 
in lArzoumanian et al.1 ( 1201 lb and iHilletalJ ( 120111) . We chose 
to do the selection manually, limiting the analysis to the main 
structures. Once these have been selected, we followed the ridge 
of the filament, connecting the column density peaks sampled 
at 20" intervals. Unless the filament has a very well defined, 
smooth ridge, this decision affects the derived filament widths. 
The structures are often fragmented and the column density 
peaks do not necessarily follow a single line. If one draws a wig- 
gly filament through every individual clump, one is measuring 
the sizes of those smaller structures. One could equally describe 
the system as a more straight filament, allowing the spatial scat- 
ter of the substructures to be counted as part of its width. Another 
critical parameter is the spatial extent analysed. For the FWHM 
values the effect is obvious and the values increase with the size 
of the analysed area, i.e., with the maximum distance from the 
centre of the filament. This is caused by the hierarchical struc- 
ture of the ISM and the need to apply background subtraction 
at some scale to remove the contribution of extended emission 
that is not directly connected with the filament. The fit of the 
Plummer profile may be less affected, but only as far as the ac- 
tual column density profile can be precisely fitted with that par- 
ticular functional form. Because of these ambiguities, the width 
of the filament and its other parameters are not well defined and 
different studies can be compared only when the methods and 
the probed linear scales are identical. In particular, to be able to 
compare sources at different distances, we were forced to select 
+0.4 pc as the extent of the structures examined. 

There are several reasons why our dust observations may 
give an incomplete picture of the structure of the clouds. Firstly, 
the column density values may be biased in the presence of tem- 
perature variations that exist along the line of sight and within 
the beam. Because warm dust emits more strongly, the colour 
temperature overestimates the true, mass weighted dust tempera- 
ture. Consequently, the column densities will be underestimated 
and the effect is most notable towards the centre of the filaments 
where the temperature variations along the line of sight are the 
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Fig. 15. The fields sorted in increasing number density of the 22 yum sources in the high column density areas relative to the average 
source density of the field. The circles are the excess of sources per square degree. The horizontal dashes show the total source 
density, including both the high and the low column density areas. 



largest. In our analysis of Sect. 13.1.41 we use column density 
maps that are based on lower, 40" resolution colour temperature 
information. To quantify the associated biases, we examined ra- 
diative transfer models of cylinders with Plummer like density 
profile (lYsard et al.l l2012h . The dust temperature distributions 
were solved assuming that the cylinders are heated externally by 
an interstellar radi ation field with intensity applicable to the solar 
neighbourhood (Mathi s et aDl98 3). The resulting synthetic sur- 
face brightness maps were analysed as the SPIRE observations, 
deriving colour temperature maps at the 40" resolution and the 
column density maps at the 20" resolution. The data were then 
fitted with Plummer functions and the obtained parameters were 
compared to the actual values of the column density distribu- 
tion. For a filament pc - 10^ cm"-', 7?gat=0.58 pc, and p = 2.0, 
the correct parameters were recovered with an accuracy better 
than 10%, even when the source was moved to a 1 kpc distance 
and a 3% observational noise was added to the surface brightness 
values. The relative errors increase with the density of the fila- 
ment, also because the density peak is less well resolved. With 
Pc = 10^ cm"-', 7?gat=0.0055pc, and p - 2.0, the parameter pc 
was underestimated by ~3Q% and Rfi^i similarly overestimated 
by ~3Q% for a source at lOOpc. For a more distant target with 
(Ikpc), the typical error was a factor of three larger. The es- 
timates of the parameter p still remain correct to within 30% 
in all the tests. The results suggest that the central density will 
be underestimated and the filament radius, TJflat, overestimated 
but, for the typical filaments in our sample, the errors are below 
50%. In the tests the density of the filaments perfectly followed 
the Plummer profile. Therefore, it does not indicate if the re- 
sults would be sensitive to small deviations from the Plummer 
shape such as, e.g., in the presence of asymmetries or substruc- 
ture. However, it is clear that the pc and Rn-^t parameters are in- 
trinsically anticorrelated. 



Further studies are required to find out how the combina- 
tion of low resolution colour temperature and higher resolution 
surface brightness data compares with the normal procedure of 
convolving all data to the lowest common resolution. 



4.3.2. Filament Structure 

The width of the filaments (or more generally, of the selected 
structures shown in Appendix ICli were found to cover a large 
range from ~0.1pc to close to 1 parsec. 

The FWHM widths are often larger than 0. 1 pc and also 
vary along the length of the filaments. The width is often an- 
ticorrelated with the column density, the smallest FWHM values 
being associated with the highest column densities, i.e., dense 
clumps and cores. Part of this anticorrelation is caused by the 
non-continuous nature of some filaments where no real column 
density excess is detected between the clumps. Depending on 
the structure of the surrounding ISM, the widths obtained in the 
gaps between the clumps can be very large. However, there is 
also real anticorrelation that can be seen in fields like G300.86- 
9.00 (PCC550), G167.20-8.69, and G227.95-2.98. 

The distances between denser cores can be as low as ~0. 1 pc 
(e.g. Gl 10.89-2.78) but are typically 0.5 pc or higher and thus 
significantly above the Jeans length. For example, G167.20-8.69 
shows a very regular pattern with column de nsity peaks a t 0.6 pc 
intervals. The results are consistent with Henning et alj (1201 Oh 
who measured a core separation of 0.9 pc in the infrared dark 
cloud filament GO 1 1 . 1 1 -0. 1 2. 

Several fields are actively forming stars and many of the ob- 
served d ense cores are likely to be g r avitatio nally bound as dis- 
cussed in lPlanck Collaboration et aTl (|201 Ibh . This strongly sug- 
gests that many of the mapped cores are undergoing gravitational 
collapse. However, anticorrelation between column density and 
filament width is expected even for isothermal gas cylinders in 
hydrostatic equilibrium. Similarly, when hydrostatic spherical 
clumps reside inside hydrostatic filaments, the former are likely 
to have smaller FWHM values. Figure [16] shows Plummer pro- 
files and the column density profiles for Bonner-Ebert spheres 
with masses equal to the filament mass within one Jeans length 
(assuming r=10K). The cores are seen to have higher column 
densities and smaller FWHM values. Note that in this calcula- 
tion, if the central density of the filament is very high, the peak 
column density of the corresponding Bonnor-Ebert sphere could 
remain lower because of the very short Jeans length. The fig- 
ure also illustrates that if the filament profiles extend far beyond 
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the 0.4 pc distance used in our analysis, the FWHM estimates 
(as opposed to the values obtained from the profile fits) will be 
biased toward lower values. 

A few sources (G89.65-7.02, G159.34+11.21, G176.27- 
2.09, G276.78+1.75) could not be fitted properly using the lim- 
ited range of ±0.4 pc. In most cases this results from the fact 
that one is looking at some elongated structures whose basic na- 
ture is different from that of actual filaments. However, the fields 
G176.27-2.09 and G276.78h-1.75 clearly do have the appearance 
of regular filaments while their widths are closer to 1 pc instead 
of the more typical ~0.2-0.3 pc. In both cases the estimated dis- 
tances are large, ~2kpc. This again raises the question whether 
these could in fact be more nearby clouds. 

The widths o f our filaments are on the average larger than 
those reported bv' Arzoumanian et al.l (1201 1 ) for the filaments in 
ICS 146. In some cases our structures might be further divided 
into smaller filaments, depending on the way th e filament de- 
tection is performed (cf. lArzoumanian et al.ll201 1', Fig. 1). Thus, 
clouds can exhibit a hierarchical structure also regarding the fil- 
aments. 



4.4. Coreshine 

Out of the 56 fields examined for the presence of coreshine, 
12 are too close to the galactic plane and show background 
emission above the threshold. They are not considered any fur- 
ther At 12yum, 9 fields showed (weak) emission (most prob- 
ably linked to PAH emission), 13 fields were not detected 
and 22 showed absorption. Emission at 3.4//m was found ei- 
ther in correlation with the emission at 12/im (PAHs) and 
in a few cases in correlation with absorption at 12yum (core- 
shine). Apart from the well-known L183 case (G6. 03 -1-36.73), 
coreshine was clearly identified in three other fields only, 
G4. 18+35.79 (LDN134), G210.90-36.55 (LDN1642, MBM20- 
21) and G300.86-9.00 (PCC550, in the Musca Complex). 
L DN 1 34 was e rroneo usly reported as showing no coreshine 
in iPagani et all (l20Tol) . LDN 1780 (G358.96+36.75), the third 
cloud in the same complex as LDN 134 and LDN 183, is bright 
at 12jt/m, unlike the two others and is therefore dominated 
by PAH emission. There are six additional fields with clumps 
showing possible signs of coreshine (G130.37H-1 1.26, G167.20- 
8.69, G173.43-5.44, G181. 84-18.46, G215.44-16.38, G298.31- 
13.05). However, the MIR intensity ratio is seen to have large 
variations that are not correlated with the dust column den- 
sity (or the dust temperature). Therefore, these detections can 
be said to be only tentative. The sample also includes targets 
like G39. 65-1-1.75 where the relative intensity of the 3.4 /im 
and 4.6 fim appears to be determined mainly by the extinc- 
tion of background radiation. Most of the fields with tentative 
or positive coreshine detection are at high latitudes. Although 
not an absolute requirement, a low background surface bright- 
ness makes it easier to detect the faint MIR excess caused by 
the scattering. In many other medium to high latitude sources 
(e.g., G108.28+16.68, G212.07-15.21, G315.88-21.44) the re- 
sults were negative or inconclusive. Finally, only 7 to 18% of the 
fields show confirmed or possible coreshine. Th i s is m arkedly 
lower than the results reported by iPagani et al.l (1201 Ol) with a 
~50% detection rate. 

It is not yet understood why some cores show the core- 
shine effect while others do not. Time is needed to let the grains 
grow but time is clearly not the only factor A high peak col- 
umn density also is not a necessary nor a sufficient criterion 
as coreshine is clearly detected towards sources with N(li2.peak) 



< 1 X 10^^ cm"^ and still not towards sources at similarly high 
latitude but with peak column densities 2-3 times as high. The 
low resolution of the WISE experiment is possibly a strong lim- 
itation, combined with the small extent of many of the clumps. 
The coreshine studies would benefit not only from higher spatial 
resolution but also from data with a higher signal-to-noise ra- 
tio. The Spitzer program Hunting Coreshines with Spitzer (PI R. 
Paladini) is currently carrying out observations that, compared 
to the WISE data, will be more sensitive by a factor of ten. The 
targets of that survey, 90 in number, have been selected from the 
Planck Early Cold Cores Catalogue (ECC) and correspond to 
~10% of the ECC. Detailed modelling is needed to separate the 
coreshine from the effects of the attenuation of the background 
radiation and to quantify the degree of the grain growth. 

5. Conclusions 

We have examined a sample of 71 fields that were mapped with 
the Herschel SPIRE instrument as part of the key programme 
Galactic Cold Cores. The examination of the sub-millimetre 
Herschel observations and the available AKARI and WISE in- 
frared data leads to the following conclusions: 

- The data confirm the presence of cold dust with colour tem- 
peratures of the total intensity typically going down to ~ 14 K 
or below. With the subtraction of the local background emis- 
sion, the estimated temperature of the larger clumps can be 
at least 1-2 K lower, the eff'ect being larger if one assumes 
that the dust emissivity spectral index in these regions rises 
above p -2.0. 

- About 50% of the fields have a mainly filamentary struc- 
ture, the filaments often being fragmented into a number of 
clumps. In several cases the morphology belies direct dy- 
namic interaction, e.g., in the form of sharp interface layers. 
Several clearly cometary clouds are also seen, again com- 
bined with a significant amount of smaller scale structures. 

- The fields include a few quiescent clouds but most clouds 
show some signs of ongoing star formation. Young stellar 
objects, as detected in mid-infrared data, are found to co- 
exist with clumps having cold sub-millimetre spectra. 

- The mid-infrared data were searched for signs of coreshine, 
enhanced mid-infrared scattering potentially caused by the 
increased grain sizes. Clear signal was observed in four fields 
with a further six tentative detections. 

- The main filamentary structures found in the fields were 
analysed. The typical filament width was observed to be 
~0.2-0.3pc. The widths were anticorrelated with the col- 
umn density, a behaviour expected both for isothermal, hy- 
drostatic filaments of different densities and for gravitation- 
ally bound spherical cores embedded in filaments. 

- Statistically, no anticorrelation was found between the mid- 
infrared emission and the large grain temperature, as ex- 
pected in the dust coagulation scenario. This can be ex- 
plained by the fact that these wavelengths probe entirely sep- 
arate parts of the clumps. On the other hand, several clumps 
were opaque enough to be seen as mid-infrared dark clouds, 
with a sharp decrease in the 22 fim surface brightness. 

Acknowledgements. MJ, JM and NY acknowledge the support of the Academy 
of Finland Grants No. 127015 and 250741. NY acknowledges the support 
of a CNES post-doctoral research grant. This publication makes use of data 
products from the Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the hifrared Processing and Analysis 
Center/California Institute of Technology, funded by the National Aeronautics 
and Space Administration and the National Science Foundation. This research is 
based in part on observations with AKARI, a JAXA project with the pailicipation 
ofESA. 



17 



M. Juvela et al.: Galactic cold cores 111. General cloud properties 



1^ 10"' 



10 



-3 



1 1. 

a. 


1 1. 

'*nM=1.0 

;- * ^ M = 3.7 




1 I- 


; ~ ^^^^S4/8.2 

"-"^ 2/0.05^-\Nl 

1 ^ I- 



10"^ 

R ipc) 



10^ 



10"^ 

R ipc) 



10' 



Fig. 16. Radial column density profiles of selected Plummer profile filaments (solid lines) and Bonnor-Ebert spheres (dashed lines). 
a) Two filaments with 7?flat = 0.1 pc and p -2 and a central density of pc - 10^ cm"^ or pc = 5 • 10^ cm"^, and a Bonner-Ebert 
sphere with a mass of 6.0Mq that in both cases corresponds to the filament mass within one Jeans length, b) Filaments with a central 
density pc -10^ cm"^ with values p/Rnm as indicated in the figure. The dashed lines are column density profiles for Bonnor-Ebert 
spheres with mass equal to the mass of the filaments within a distance of one Jeans length. 
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Table 1. The observed fields. The columns give the field name, the centre po- 
sition in Galactic and equatorial coordinates, the size of the mapped area, the 
distance, flags for the source morphology, and a list of possible other identifica- 
tions. 
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^ The r eference for the distance estimate is given in parentheses: 1 - association with 'Simon et al.' (2006) IRD C, 2 - reddening of SDSS stars 
(seelMc Gehee 2012 in preparationh ). 3 - 3D extinction mapping using 2MASS catalogue stars (see Marshall et al.ll200^ . 4 - kinematic distance, 
5 - association with a known molecular cloud complex, 6 - direct identification with an object with a distance estimate available in the literature. 

* Characterisation of the field morphology: F - filamentary, B - boundary, C ~ cometary, 1 - isolated, X - complex. 
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Table 1. continued. 
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Table 2. Masses of the mapped fields and of the selected clumps. The columns are 
(1) the name of the field, (2) the total mass of the mapped area, (3)-(5) coordinates 
and radius of a selected aperture, (6)-(8) temperature, peak column density, and 
mass within the aperture as derived from the total surface brightness without 
the subtraction of the local background, (9)-(10) temperature and mass based 
aperture photometry with a local background subtraction. 
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M. Juvela et al.: Galactic cold cores III. General cloud properties 



Table 3. Properties of selected filaments. The columns are: (1) name of the field, 
(2) the length of the analysed filament, (3) width of the analysed filament area, 
(4) total mass of this area, (5)-(6) average column density along the filament and 
its dispersion, (7) Jeans length calculated for average and peak column density, 
(8) measured FWHM of the mean profile of the filament, (9)-(ll) parameters 
of the Plummer profile fit, (12) mass per unit length according to the fit of the 
Plummer function. 
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Appendix A: Maps of observed and derived quantities 

Figures ??-?? show maps of observed and derived quantities for the fields other than the one already shown in Figs.|4]-[5] Included in 
the figures are the colour temperature map derived from SPIRE data (frame a, 40" resolution), the 250 SPIRE surface brightness 
map (frame b, 18" resolution), the AKARI wide filter maps at 140 jum and 90yum (frames c and e; convolved to a resolution of 
one arc minute), the visual extinction Ay derived from 2MASS catalogue stars (frame d; resolution of 2'), and the WISE 22 //m 
intensity (frame /; 12" resolution). In the cases where public WISE data were not yet available, frame / shows the 25 jjm IRAS 
map (resolution of 4.5')- The respective beam sizes are indicated in the lower right hand corner of each frame. The black circles 
indicate the positions of clumps selected for further examination in Sect. 13.2.11 and the white circles indicates a reference region 
used for background subtraction in the SED plots in Sect. l3.2.TI In frame a, the arrow indicates the position of the stripe plotted in 
AppendixlB] 

Appendix B: Column density maps 

The figure ??-?? shows the column density maps for all the fields excluding those already shown in Fig. |6] 

Appendix B: One-dimensional cuts of the surface brightness data 

Figures ??-?? show selected surface brightness values along the one-dimensional cuts marked in the frames a of Fig.|4]and Figs. ??- 
??. 

Appendix C: Figures of the selected elongated cloud structures 

The analysis of the filamentS in the fields G163.82-8.44 and G300.86-9.00 (PCC550) were shown in Figs.EHE] The plots for the 
other fields listed in Table[3]are shown in Figs. ??-??. These include only the fields with an existing distance estimate and where a 
clear filament or other distinct elongated structure could be discerned. 

Appendix D: Spectra of selected regions 

The spectral energy distributions of selected regions were presented in Fig.[TOl Similar figures for the remaining fields are shown in 
Figs. ??-??. 



